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Given two spherically symmetric and short range potentials Vq and Vi 
, for which the radial Schrodinger equation can be solved explicitely at zero 

energy, we show how to construct a new potential V for which the radial 
equation can again be solved explicitely at zero energy. The new potential 
and its corresponding wave function are given explicitely in terms of Vq and 
\ Vi, and their corresponding wave functions (po and tpi. Vq must be such that it 

' sustains no bound states (either repulsive, or attractive but weak). However, 

Vi can sustain any (finite) number of bound states. The new potential V 
has the same number of bound states, by construction, but the corresponding 
(negative) energies are, of course, different. Once this is achieved, one can 
start then from Vq and V , and construct a new potential V for which the 
5^ \ radial equation is again solvable explicitely. And the process can be repeated 

indefinitely. We exhibit first the construction, and the proof of its validity, 
for regular short range potentials, i.e. those for which rVo(r) and rVi(r) are 
at the origin. It is then seen that the construction extends automatically 
to potentials which are singular at r = 0. It can also be extended to Vq long 
range (Coulomb, etc.). We give finally several explicit examples. 
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I. Introduction 



Consider the reduced radial Schrodinger equation for a spherically symmetric 
potential V{r) [T] 
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V{r) 



r G [0,oo) , k>0 , (fie{k,0) = , 
V{r) real , locally for r 7^ , and l^(oo) = . 



(1) 



We exclude therefore confining potentials like the harmonic oscillator, etc. 

There are only a few potentials for which the radial Schrodinger equation can 
be solved explicitely for all k and all i. These are essentially the square well, the 
Coulomb potential, sums of 5-function potentials pQ ; and for potentials which are 
more singular than at the origin, but repulsive, only A/r^ [2^, for which the 
solution is given in terms of very complicated Mathieu function. 

If we restrict ourselves to the case of one single i, then we can include the 
Bargmann potentials ^1 E] , for which, by construction, the radial Schrodinger equa- 
tion can be solved for one specific value of i, and for all k. We remind the reader 
that Bargmann potentials are those for which the S-matrix Se{k) is a meromorphic 
function of k in the fc-plane. They can be constructed for each i. 

In the particular case of £ = 0, the radial equation can be solved for all k for the 
following potentials P : 



VAr) = Xe-^' 



1/2 (r) = — (Hulthen) 



1 - e-^'^' 



(2) 
(3) 



and, more generally, the Eckart potentials |1] : 



T/ ( \ ^1 e-^'- ^ A2 e-^^ 



1 + C e-^^ (1 + C e-'''- 



i2 ' 



(4) 
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of which (12)) is a particular case. The solutions are given in terms of hypergeometric 
functions. 

In the case of = 0, and for all one can add the potential [3] : 



VAr) 



A r 



a-2 



, a > , ^ > 



(5) 



An interesting particular case is when a = 2 



A 



2 ' 



> . 



(6) 



Finally, for /c = and £ = 0, one can solve the radial equation also for 



VAr) 



(C + r) 



. 



(7) 



We shall give later the explicit solutions for some of these potentials, when they are 
simple. 

The purpose of the present paper is to show that if one can solve explicitely the 
radial Schrodinger equation at = for 



Voir) + 
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J r|Vo(r) dr < oo , j r^^^'^\Vo{r)\dr < oo 



no bound states, 



and for k = and £ = for the potential 



r\Vi{r)dr < oo , / r \Vi{r)\dr < oo 
Ji 



any number of bound states. 



(9) 



then one can solve explicitely the radial equation, again at = and £ = 0, for the 
potential 



V{r) = Vo{r) + Fo{r)Vi{Go{r)) , (10) 

in terms of the solutions for (jH)) and ©• Here, the functions -Fo(r) and Go{r) are 
given very simply in terms of the solutions for (jSI). 

To begin with, we consider the case i = 0. If we call by (po{r) and Xoi^") the two 
independent solutions for Vo(r) defined by 

r ^o(O) = , <f',{0) = 1 ; xo(0) = 1 , 

(11) 

Fo(r) and Go{r) are given by 

Fo{r) = [xoir]]-' , Go{r) = ^ . (12) 

As we shall see later, since by assumption, Vo(r) has no bound states, all the above 
quantities are meaningful because both (po{r), and Xo{r) defined by 



Xoir) = ipo{r) I , (13) 



do not vanish anywhere, except for ipo ai r = (see details in the next section). 
Note that the converse of (fT^ is 

rr 

Mr) = Xoir) ^dt. (13') 

The solution of the radial equation at = and i = for the potential V{r), 
ffTUI) . is then given by 
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(^(r) = xo{r) <f\ 



[ ifiO) = , 



(14) 



where ipi is the regular solution of the radial equation for Vi, Q, defined by 
(pi{0) = 0. By assumption, both ipo and (pi are known, and both vanish at r = by 
definition. The above formula can be checked directly by differentiation. We shall 
see in the next section how it was found. 

Remark 1. As we shall see, x = maps r G [0, oo) into x G [0, oo). The 

mapping is one to one because of below, and is, of course, twice differentiable 
(see dH))). 

Once (fT^ is known, it is easy to generalize it to the case one has angular moment 
with Vb(r) : 



Voir) + 



> 



rVoir) G L\0) and r^^+^Vo{r) G ^^(oo) , ^ 
no bound states, 

Vi{r) remaining unchanged, (fo and Xo are now the solutions of the radial equation 
with the potentials (8'), so that (fTUI) becomes 



V{r) 



Voir) + 



+ 1) 



+ Fo(r)Vi(Go(r)) . 



All other formula given above remain unchanged. 

In short, if one can solve explicitely the radial equations 



(100 



Voir) 



Voir) + 



iii+l] 



'Voir) , V5o(0) = 



^';ir) = V,ir) y.i(r), ^i(O) = 



(15) 
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where Vq and Vi satisfy the assumptions shown in (8') and Q, then the solution of 



with ip{0) = 0, is given by 



XW) \Xo{r)J 



^(r) = Xo(r) ^1 f ^) . (17) 

Remember that Xo('") is always defined by (13). 

It is easy to check our assertion by differentiating twice given by p7|) . and 
using (fT3j) . We shall see in the next section how (fl^j) and (^7j) were found. We shall 
also see that one can replace Vo(r) by strongly repulsive potentials which are more 
singular than at the origin. Examples will be provided for 

Vo{r) = ^, g>0, n>2, (18) 
for which the radial Schrodinger equation is soluble for all £ at = [H]. 

Remark 2, the bound states. As is well-known, the nodal theorem [Zj asserts 
that the number of bound states of V{r), ()10|) or (10'), is given by the number of 
the nodes of the regular wave function (p{r), ()17|) . Since neither Xo{^) for ^ ^ 0, nor 
(po{r) for r > 0, do not vanish (remember that, by assumption, Vq has no bound 
states), and x = ^^44 maps r G [0, oo) into x G [0, oo) and the mapping is one to 
one, it is obvious on (fTTjl that and ipi have the same number of nodes. Therefore, 
Vi, (P), and V, (fTn|) or (10'), have the same number of bound states. Of course, 
the energies of these states are different for Vi and V. In any case, one has also the 
Bargmann bound for the number of bound states PI El El ^ 



oo 



n{V)=n{Vi)< r\Vi{r)\dr < oo . (19) 







Remark 3, iterating the process. Once we have the explicit solution ()17|) 
for the equation (fTBj) . we can start now with the couple [Vo(r), V"(r)], instead of 
[Vo(r), Vi(r)], and look for the solution of the radial equation at A; = for 
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V{r) 



Voir) + 



+ 1) 



+ 



■V 



Xo{r) 



We will find now, of course, the solution 



And this process can be continued as many times as we wish. 



(20) 



(21) 



We end this introduction by giving one example with the potentials 



Voir) 



Mr) 



;i + ar)4 ' 
/ 1 + ar 



Xoir) = (1 + ar) 



sinh 



cosh 



a > , A > 
^ \/X ar\ 



^ 1 + ar ^ 
-\/A ar\ 



cosh . , / \/X ar ' 
smh 



1 + ar J sinh a/A \1 + ar ^ 



and 



Viir) 



Mr) 



(62 + r2)i/2 



=— sinh ( V5' - 1 Arctg - j 



(22) 



(23) 



from which one can calculate by formula (|T7jl . Now, since the Schrodinger equation 
can be solved for the potential (j23|l for all £, we can invert the roles of Vq and Vi, 
and start with 



VM\ 
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^2 I ^2^ 



(24) 



For general £, the solutions <y9o and xo are given in terms of hypergeometric functions 
jSj. We restrict ourselves to the case of £ = 0, for which 
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y/W+r^ . , / / . r 

^y^— y smh l^^g - 1 Arctg - 



Xoir) = A/P~+T^cosh ( y/g — 1 Arctg - 



(25) 



and take (|22|1 as the second potential, with 



1 + ar ( \fX ar' 
ipi{r) = ^ sinh 



1 + ar 



(26) 



II. Derivation of the solution (jl4p 

We begin by studying the properties of the solutions y:>o{r), and Xo(^) defined by 
(Hn)), of the radial equation at k = and i = : 



' foi'^) = Mr) ^o{r) , 
. <^o(0) = , xo(0) = 1 



(27) 



where Vo(r) satisfies the assumptions shown in (jH)). Since the radial equation is 
homogeneous, we can normalize its solution ip as we wish. For ()27j) . the usual 
convention is to put 



<^o(0) = 1 • 



(28) 



Then we can combine ()27|1 and ()28j) into a single Volterra integral equation PI El E] 



(po{r) = r + (r - r')Vo{r')(po{r')dr' 
Jo 



(29) 



It can then be shown that, iterating the above equation, and using the assumptions 
on Vo(r), namely that rV is at r = 0, r^V"(r) is at r = oo, one gets an 
absolutely and uniformly convergent series defining the solution (po{r), together with 
the bound P H 
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(r)| < r e/o'''-'l^"('-')l'^^' < Cr , (30) 



where C is an absolute constant less than exp /q°° r|Vo(r)|c?r. Using this bound in 
()29|1 . we find that indeed, for r — > 0, we have 

¥Po(0) = , ¥''o(0) = l, (31) 

and for r oo, 



^o{r) = r 







r'Vo{r')ipoir')dr' + o(l) . (32) 



1 + / Vo{r')Mr')dr' - 
Jo 

where all integrals are absolutely convergent. 
There are now two cases : 

i) the potential Vo(r) is positive. Then it is obvious on the iterated series of 
that all the terms are positive, and so is v^o(^) for K follows then from ()27|) 

that ipo{r) is a positive convex function of r. It increases indefinitely, and we have, 
on the basis of 



(po{r) > , and convex, 

(33) 

ipo{r) = Ar + B + o(l) , asr^oo, l<yl<oo, B <0 . 

A schematic picture of cpo is shown on Fig. 1. 

ii) Vo(r) < 0, but not strong enough to have bound states. Then we find from 
the nodal theorem relating the bound states of Vo(r) to the zeros (nodes) of ipo{r) 
for r > [7|, that (po{r) is again positive, and since Vq is now negative, (po is positive 
and concave. A schematic picture of ipo is shown on Fig. 1. From (jHSI), is now 
replaced by 

(po{r) > , and concave, 

(34) 

y9o(r) = Ar + B + o(l) , asr^oo, 0<v4<l, B >0 . 

Remark 4. Here, if A < 0, then, since ^Pq{0) = 1, and (y9Q(oo) < 0, ipo must have 
a zero in between, and therefore one has a bound state, in contradiction with our 
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assumption of no bound states. If A = 0, this means that one is at the threshold 
of having a bound state. More precisely, that one has a resonance at zero energy 
in El El 5 a possibility we have excluded also. 



We come now to the second, independent solution Xo{r), defined by (fTT?j) . First 
of all, since fo{^) is always positive for r > 0, and from and ((MI), the integral is 
absolutely convergent at its upper limit, and so Xo is twice differentiable, and satisfies 
the same equation as (fQ. It is trivial to show that the Wronskian of the two, (lllj). 
is 1. Now, when r — > 0, the integral in (jl3p diverges, but since (po{r) = r + o(l) as 
r — i> 0, and there is ipo{r) in front of the integral, it is trivial to show that we have 
Xo(0) = 1, as shown in Also, on the basis of the Wronskian (jlip . and (j3ip . we 



This general property is a consequence of rVo(r) G at the origin. The derivative 
of Xoi'f') at r = maybe finite or infinite, depending on the behavior of Vo(r) near 
r = 0. If Vo(r) itself is at r = 0, one can write also the integral equation IHEIIII : 



and iterate it, as we did with ()29|) for ipo, to find the solution, which turns out now 
to be bounded everywhere. One then immediately sees on (IHUj) that Xo(0) is finite. 
We have, therefore, according to ^T^. and or dSl (see Fig. 2) : 

' Xo(0) = 1 , Xo{r) > for all r , 
Xo{^) is a convex and decreasing function when Vq > , 



find 



lini r Xo{r) = . 
1 — >o 



(35) 




(36) 



Xo{r) is a concave and increasing function when Vq < 0, 
with no bound states. 



(37) 



Xo(oo) = -^0,^;Vo>O^A>l,Vo<O^A<l. 
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Consider now the mapping 



x{r) 



Xo{r) 



(38) 



According to (j37|) . this a perfectly regular and differentiable mapping, and is one to 
one since, according to (jllj) . we have 



dr Xoir) 



Xoir) 



> . 



(39) 



It follows then, since (poir oo) — > oo, and xoir oo) ^ ^ 7^ 0, oo, that the 
mapping is one to one : 



r G [0, 00) -v^ X G [0, 00) , x(0) = , x(oo) = 00 
In fact, this mapping is twice continuously differentiable since 



(40) 



-2Xo(r) 
Xoir) 



(41) 



and Xo('") is a continuous function of r for r > 0. This last property follows from 
x'oir) = yoir)xoir), where, by assumption, Vo(r) G for r > 0. Since Xo('^) is a 
continuous function, Xo('^) is also for all r > 0, and so Xo(^) is continuous for 
r > 0. Xo cannot have jumps 0. See Fig. 3. 

Once we have established that the mapping (^01) is regular and twice continuously 
differentiable, we can consider the equation 



^"ir) = {V^{r) + Viir)]^{r) 
y.(0) = , (/.'(O) = 1 , 



(42) 



where Vq and Vi satisfy the conditions shown in (jHJ and Q. We make now the 
change of variable and function 



X 



Xoir) ' 



Tpix) 



'^ir) 
Xoir) 



(43) 



r=r{x) 
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where r{x) is the inverse mapping, i.e. the inverse function of x = x{r). Obvi- 
ously, r{x) is also twice continuously different iable. Differentiating now twice ip 
with respect to x, and using (fT^ . we easily find 



'4'{x) = Xoir)Vi{r) _ i;{x) 



r=r{x) 



(44) 



There is no longer Vq present. From the definition of of if{r) given in (|42j) . 

and (jnSl), it is obvious that, because (p{r) = r + o(l) as r — > 0, we have 



iPiO) = , ^(0) = \imjj{x) = lim [<^'(r)xo(r) - <^(r)x[,(r)] = 1 . (45) 



Suppose now that, from the beginning, Vi(r) in (j42p was of the form. 



1 vJ,-^»(^) 



Xo{r 



Xo[r) 



(46) 



where xVi{x) E at x = 0, and x^Vi(x) E aX x = oo. Then (jl^ would become 



iIj{x) = Vi{x) 'ip{x) , 

^ V^(o) = , ^(0) = 



(47) 



which we assume to be explicitely solvable. Then, from the definition (j43|) . we would 
have for the solution of (j42p . with a Vi of the form (j46|) . 



Mr) 



(48) 



Combining all these, with a slightly different notation, we have therefore the follow- 
ing 



Theorem 1. The solution of 
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¥^(0) = , ^'{0) = 1 , 



where ipo and xo are the two solutions of 



' foi^) = Mr) Mr) , 

V3o(0) = , v^o(O) = 1 , no bound states, 
[ Xo(r) defined by (HI, xo(0) = l, 



is given by 



V{r) = xo{r) 



Mr) 
.Xo{r) 



where (pi{x) is the (regular) solution of 



' M^) = Vii^) 'fii^) ) 
^ <^i(0) = 0, 0i(O) = l. 



(49) 



(50) 



(51) 



(52) 



Therefore, if the Schrodinger equation at k = and i = can be explicitely solved 
for Vo and Vi, then the solution of ()49|) is of the form (j^T|) . As we said in the intro- 
duction, one can check directly, that (jM]) is indeed the solution of (jl^ . For bound 
states in and (j^ . see below, after Remark 6. 



Higher waves, i > 0. So far, we have been assuming i = 0. It is easy to extend 
the results to the case £ > in (|^. i.e. to begin with 



Mr) 



Voir) + 



iii + 1) 



Mr) 



(53) 



and then add the potential Vi in (pUj) . Vq is, as before, supposed to be such that 
tVq G at r = 0, and r^Vo G at r = oo. We shall see later that we need more 
rapid decrease at infinity. The regular solution ipQ is usually normalized as follows 
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One can then combine (jH!^ and (jH^ into the single Volterra integral equation 

Solving this equation by iteration, one finds again, as for the case £ = 0, an absolutely 
and uniformly convergent series defining the solution (pQ, together with a bound 
similar to (j3(jp for all finite r > : 

\Mr)\ < C /+^exp (^J\'\Vo{r)\dr'^ < C , (56) 

where C and C are absolute finite constants [HElIll- Using in (|^. one sees 
immediately that, for r — > 0, 



For all the above results, we need only rVo G -^^^(0). Also, by assumption, there 
are no bound states for (j^Hj) . It follows again that, by the nodal theorem (7j, ^^ir) 
cannot vanish for r > 0. Because of ()57|). we find therefore that 

(/?o(r) > for all r > . (58) 

From this, and ()53p. it follows immediately that if Vo(r) > 0, then V5o(^) is convex. 
For Vo(r) < 0, the situation is more subtle than for the case of £ = 0, nd y^o may 
become concave in some interval [Rx^R-i). 

We wish now to look at the behaviour of v5o(^) as r ^ oo. We assume here 

r''+' V^{r) e L^(oo) . (59) 
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Then we can let r ^ cxd in ()55|) . to find 

1 1 

9^0 



(2£+l)!! Jo 2£ + l 



r'-%{r')^Q{r')dr' 



2i + l 



(60) 



Since '{)o{r) never vanishes (absence of bound states), must be always positive. 
If Vo > 0, then Ai > 1, and if Vq < 0, < < 1. For B^, it is just the opposite. 
These are quite the same as for i = 0. 

We can now define the second (independent) solution Xo(^) by (fr5|) again, and 
we find now, using (jHTjl and (jHUjl . that 



Xo[r) 



Xo[r) 



(2^-1)!! 



+ 



r ^ 



(61) 



i2£+l)A 



-r , r ^ oo 



If we introduce the same variable x = x{r) as before 



X = x{r) 



Xo{r) 



(62) 



we find 



x{r) = A r2^+^ H 

x{r) = B r'^^^^ J r 



r ^ 



(63) 



oo . 



The rest of the analysis goes as before, and we find 



Theorem 2. Theorem 1 is vahd if we add i{i + l)/r^ to Vq in (jl9|) and (j50|) . i.e. 



/'(r) = 
(^(r ^ 0) 



K)(r) 



(2£+l)!! 



+ ■ ■ 



(64) 



and 
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/oW= K,(r) + 



+1) 



(/?o(r) , r 



2^+2 Vo(r) G Li 



(oo) , 



< 



(65) 



<^o(r 




+ ■ ■ ■ , no bound states, 



(2£+ 1)!! 



Xoir) defined by xoi^ 



0) 



(2^-1)!! 



while remains unchanged. The 



r 



solution is again provided by ()51|) . This can also be checked directly. 

Remark 5. Since the behaviour of ^pi{x) is x+- ■ ■ as x 0, and ^Po{^)/Xo{^) = ^ 
as r — > 0, it is obvious on that we have, as we should. 



Bound States. So far, we have assumed that there are no bound states in (j52j) . 
If there are n bound states with Vi{x), i.e. in (j52j) . which is the same for Theorem 
1 and Theorem 2, this means that (fi has n nodes for x > 0. And we have also n 
nodes for the full solution (f, given always by (|5ip. The potentials Vi and V have 
the same number of bound states. But, of course, the binding energies are different. 

Remark 6. As we said in Remark 3, once we have V{r) and v^(r), we can now 
start again with Vq and V instead of Vq and Vi, and proceed as before. This process 
can be repeated as many times as we wish, and we get more and more potentials 
for which the radial Schrodinger equation at zero energy can be solved. Also, the 
number of bound states, if any, remains the same. Unfortunately, the potentials and 
the wave functions become quickly very complicated. However, one may ask what 
will happen at the limit of infinite repetitions. This question is certainly not easy 
to answer. 



ip{r) = ai r H , r ^ . 



(66) 



Likewise, it is easy to find 



(67) 
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Singular Potentials. As we said in tlie introduction, the radial equation can be 
solved at A; = and for all C for singular potentials which are just inverse powers 
potentials shown in (fTS|) . The solutions are given in terms of modified Bessel and 
Hankel functions. We shall see explicit examples in the next section. One more class 
is given by [HI HHj 

^o(r) = -jj- ^ + — -2 + -—^ , (68) 

with p < 2, gi > 0. These potentials must be cut at r = i?o < 1 in order to avoid 
none integrable singularities at r = 1. The solution is given in terms of Whittaker 
functions jniE]. Once we know the solution ifo{f), we can proceed as before, and 
add a regular potential Vi with any (finite) number of bound states. 

We should mention here that, contrary to the case of regular potentials at the 
origin, i.e. those for which rVo(r) G at r = 0, here, because of strong singularities 
at r = 0, we find [H E] ¥?o(0) = <fo{0) = ■ ■ ■ (^{,"^(0) = ■ • ■ = 0. All the derivatives of 
(pQ vanish at r = 0. The normalization is therefore arbitrary, and cannot be made 
at r = 0. Once this is chosen (usually by the behaviour of at r = oo), then Xo(^) 
is given again by (fT^ . and we have the Wronskian W{(po, Xo) = 1- Usually, in such 
a case, it is customary to normalize (po at infinity, according to 

(po{r) = r + C + o(l) , r ^ oo , (69a) 

which entails also 

Xo{r) = r + C' + o(l) , r ^ oo . (69b) 

As we shall see on explicit examples in the next section, our procedure for generating 
new potentials can go on without modifications. 

III. Examples 

1. Regular Potentials. We have already given, in the introduction, as examples 
for the applications of Theorems 1 and 2, the solutions of the Schrodinger equation 
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for i = 0,Vo = (7), Vi = (6), or Vq = (6) together with £{£ + l)/r^, and Vi = (7). 
They are given by formulae (j^^ - ()26p . More examples are obtained by combining 
any two potentials among those given by Q to (|7j). We need only the solutions cp 
and X, the latter defined by (fT^. for these potentials. 

a. Exponential potential 

' V{r) = Ae"'^'^ ,A>0,/i>0,£ = 

V / V / (70) 

Iq and Kq are modified Bessel and Hankel functions of order zero, and the constants 
a and /5 are determined to have ^^(O) = and ^^'(0) = 1. It is then easy to show 
that, according to our general analysis of section II, we have 

^^r^oc)=Io(^^y + --- (71) 
The presence of r in (71) is due to the presence of log 2; in Kq{z) where z ^ jS]. 

b. The potential |3j) for a > 0, i > 0. The solutions ip and x are given by 
combinations of hypergeometric functions of appropriate arguments. We refer the 
reader to [S] for details. 

c. Hulthen Potential (?>), I = 0. Here also, the solutions Lp and x ^"^^ given in 
terms of appropriate hypergeometric functions F . See ^ H] for details. 

2. Singular Potentials [6j. Here, we consider only three cases. 

d. Inverse Power Potentials, i >0 : 
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V{r) 



, g>0 , n>2i + 3 



2i+l 
n-2 



2v^ 



a/2 



{n — 2)r 2 
2^ 



n-2 I ; 



(72) 



aK'2e+i 

"-2 V(^ - 2)rV 



+ 0l2i+l 
n-2 



2v^ 



— 2)r "2^ 



where Ii, and /^jy are modified Bessel and Hankel functions. We must clioose n > 
2£ + 3 in order to comply witli ()fi()|) [0] . 

Tlie parameters a and (3 are determined for having x(r) to comply with the 
asymptotic behaviours deduced from (jl3p . for r ^ and r ^ 00. One has to 
remember the Wronkian 



W 



a 



(73) 



The case £ = 0, n = 4 is particulary simple. One finds 



V{r) 
X{r) 



0, 



r— >oo * 



^sinhf^l = 1 + 



(74) 



e. Logarithmic Potentials jHl . We consider here the simplest case of ()68|) with 
p = 1, and any angular momentum £ > : 



\/(r) 



a; 



+ 



9 



Log ( - ) Log^ ' - 



9{R-r) , a > , i? < 1 , 



y9(r) = r"^/^ 



rf-2z/) 



r(i-z/-/i 
r(2z/) 



M^,,(x) 



(75) 



^ ;^(r) = [a M^,,(a;) + /5M^,_,(a;)] 
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where M^ ,^ are Whittaker functions jH], 

x = (2£ + l)Log^ , k = ^^, u = z^^^, (76) 

and a and (3 are determined according to (|T3|l for r ^ and r — » oo. Note here 
that, for r > i?, we have the free equation (no potential), and, therefore, we must 
first adjust the free solution to the interior solution at r = i?, as usual. 

Note here that the singular part of the potential is just the first potential, and 
that is why we must choose a > 0. The second potential is regular since it satisfies 
rV{r) E L^{0). We can therefore choose g ^ 0. There are several more examples of 
singular potentials for which the radial equation can be solved explicitely. We refer 
the reader for details to jH]. 

/. Coulomb Potential. The Coulomb potential is regular at the origin, and so we 
have the usual solutions ip and as defined previously. We choose, of course, the 
repulsive case. The solutions can be read off from (72) adapted to n < 2, or else, be 
obtained in the standard way PQ Ej. One has, for i = : 

( a 

V = - , a>0, £ = 0, 
r 

' y,(r) = ,/^/i(2v^) , (72') 
V a 

where Ji and Ki are the modified Bessel and Hankel functions, and similar formulae 
for i > 0. As we see here, the long range tail of the potential leads to the exponential 
growth of at r ^ oo, ~ r^^'^ exp{2y/ar) , and the exponential decrease of x(^) ~ 
r^'''^ exp(—2^/ar), to zero. This does not affect the validity of the change of variable 
r X = ip/x, etc, of section II, except that now x grows exponentially as r ^ oo. 
Note that Ji and Ki do not vanish, Ji for r > 0, and Ki for r < cx3 is 
again an increasing convex function, and x ^ decreasing convex function. The only 
difference with the short-range potentials is that, now, grows exponentially, 

so that, in (jl^ and if Vq is chosen to be the Coulomb potential, the second 
potential may seem to become infinite as r — oo (x — > oo). However, we have 



20 



always assumed Vi(r) to be short range, i.e. decreasing fast enough at infinity. It 
follows that Xo^(^)^i(3^) is again short range. For instance, if Vi(r) ~ r~^, then 
Xo''Vi(r) ~ l/(a;^Log^x), etc, i.e. xxo'^Vi is in x at x = oo. 

Other long range potentials of the form g/r"-, n < 2, can be dealt with in the 
same way by adapting (72) to n < 2, and one reaches similar conclusions as for the 
Coulomb potential. As for confining potentials like the harmonic oscillator, etc, we 
shall consider them in a separate paper. 

Concluding remarks. 

So far, all the potentials for which the radial Schrodinger equation has been 
shown to be soluble analytically in closed form have their solutions given by various 
hypergeometric functions in appropriate variables El HI E]- In fact, in many 
instances, as we have seen on examples, the hypergeometric functions simplify to 
Bessel and Hankel functions of real or imaginary arguments. The only exceptions 
are Bargmann potentials PlEl, for which the solutions are given in terms of rational 
functions of sinha^r, cosha^r, sin jSjV, and cos jSjV, j = 1, ■ ■ ■ n, where aj and jSj are 
given by the positions of the poles and zeros of the S'-matrix, and the potential 
[2], for which the solution is given in terms of Mathieu functions. 

In the present paper, as seen on and the potentials themselves have 
their arguments given by ratios of hypergeometric functions, and the solutions are 
then hypergeometric functions of ratios of hypergeometric functions, as seen on (j^Tjl . 
And this process can be repeated indefinitely, as we saw before. One may then ask 
what the potentials and their wave functions become in the limit. Also, we saw that, 
for both Vo and Vi, we can take potentials which are very singular but repulsive at 
the origin, like gr~'"' and Ar"™, g > 0, X > 0, m,n > 2. All kinds of combinations 
are therefore possible for Vq and Vi. 

It is trivial to construct infinitely many potentials for which the Schrodinger 
equation could be solved at zero energy. One can choose any positive, convex, and 
twice continously differentiable function, which is decreasing, and such that 
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Xo(0) = 1 , Xo(oo) = - , 1< A < oo 



(77) 



and write 



K)(r) 



Xo{r) 



(7J 



If Xo(^) is decreasing fast enough to A ^ at infinity, then Vo(r) is short-range, foir) 
is defined here by (13'). Example : 



(1 + a) 



a 



, a > , /3>0, n> 1 . 



(79) 



According to (78), we have 



VJr) 



anin + 



(l + a)(l + /?r)"+2xo(r) 



r— +00 



-n-2 



^0) 



However, all this is valid for this particular Vo(r). If we try to introduce a coupling 
constant A in front of Vb(r), i.e. try to solve 



X'ir) = XVo{r) x{r) , 

we usually do not find explicit solutions. This is indeed the case here. 
Next, consider 



(81) 



Xo[r) 



1 + e- 



■fir 



/i > 



(82) 



where A = 2. This leads to 



Vn(r) 



(83) 



• ■ Xo{r) 1 + e-f^- 
(po{r) is then easily calculated from (13'). In this case, one can solve (81) for any 
A since (83) is a particular case of (jH), and the solutions are given, in general, in 
terms of hypergeometric functions jl]. Only for A = 1, they reduce to the simple 
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form (82) for Xo, and the corresponding (po{r). 
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